function u = initial_data(N,x_max)
h = x_max/N;
u = [];
for i = 0:N-1
    u = [u,initial_function(i*h)];
end
u = u';